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The two-dimensional dense 0(n) loop model for n = 1 is equivalent to the bond 
percolation and for n = to the dense polymers or spanning trees. We consider 
the boundary correlations on the half space and calculate the probability Ph that a 
cluster of bonds has a single common point with the boundary. In the limit n — )• 0, 
we find an analytical expression for Pb using the generalized Kirchhoff theorem. The 
obtained exact value is surprisingly close to the value of Pb conjectured for the dense 
0(1) loop model. 
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I. INTRODUCTION 

The dense 0{n) loop model [1] is defined by drawing two lines in each elementary cell 
of the square lattice. Two possible states of the cell are shown in Fig.l. The lines on the 
whole lattice with appropriate boundary conditions form a system of closed loops with the 
Boltzmann weight n ascribed to every loop. 




FIG. 1: Elementary cells. 



A rise of interest to the 0{n) model, especially to the case n = 1, was triggered by 
remarkable observations concerning the ground state of the antiferromagnetic XXZ quantum 
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FIG. 2: Loops on the horizontally periodic lattice. 

chain and, equivalently, to properties of the transfer matrix of dense loop model taken in 
the Hamiltonian limit of extreme spatial anisotropy [2-7] . It was noted that the elements 
of the ground state of the loop Hamiltonian are equal to the cardinahty of different subsets 
of configurations of the Fully Packed Loop (FPL) model which in their turn are connected 
with classes of Alternating Sign Matrices (ASM). 

The states of the 0[n) model can be defined in terms of the connectivity condition for 
points of intersection between loops and a horizontal line cutting the loops at these points. 
Two points are connected by a link if there exists a line between them via the half space 
above the cut. For instance, imposing periodic boundary conditions in horizontal direction 
for the lattice in Fig. 2 and specifying the boundary conditions at the upper edge, we obtain 
three minimal links between points 1 and 6, 2 and 3, 4 and 5 which are points of intersection 
between loops and the bottom line of lattice belonging to upper half space. A typical 
configuration of links for a larger lattice is given in Fig.3. 

Several interesting conjectures have been obtained about the probability distribution of 
configurations of links of the dense 0(1) loop model. To be closer to the aim of this article, 
we mention just one conjecture about the number of configurations containing the minimal 
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FIG. 3: A link configuration. 



link at a given position [7] 
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for even periodic system and 
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for odd periodic system. Here Aht{N) is the number oi N x N half turn invariant ASM [8]. 

Since the total number of periodic link configurations is Aht{N), the probabihty of 
minimal link Prob{. . . ^ . . .) tends to 3/8 in the limit N — > oo. Due to the simple bijection 
between the loop configurations on the square lattice and the bond configurations on its 
sublattices, the dense 0(1) model is equivalent to to the bond percolation. Then, the 
constant 3/8 has a simple interpretation as the probability of a certain class of percolation 
clusters (see Section II). Despite its simplicity, an analytical derivation of the probability of 
minimal links as well as of other link configurations remains an open problem. 

In this work, we consider the problem of evaluation of Prob{. . . ^ . . .) in the dense 0{n) 
model for the case n — 0. In the limit n — > 0, the bulk loops disappear and the system 
of lines is converted into the model of dense polymers [9]. The bijection giving the bond 
percolation for n = 1 on a sublattice, leads in the limit n — )> to the spanning trees on the 
same sublattice. A lot is known about correlation functions of the spanning tree model due 
to the Kirchhoff theorem [10]. However, we will see in subsequent sections that evaluation 
of Prob{. . . . . .) involves non-local diagrams which need a generalization of the Kirchhoff 

determinant formula. Similar correlations appear in the Abelian sandpile theory [11-15], in 
the problem of loop-erazed random walk [16, 17], in the problem of "watermelon" embedded 
into the spanning tree [18, 19]. Below, we use the generalized determinant formula to find 
Prob{. . . . . .) and compare the result obtained for n = with that conjectured for n = 1. 
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II. LINK STRUCTURE AND BOUNDARY CLUSTERS 

The square lattice of sites with integer coordinates {n, m) can be divided into two sub- 
lattices, black and white. For sites of the black sublattice n + m is even, for sites of the 
white one n + m is odd. The bijection between loop and bond configurations mentioned in 
Introduction is shown in Fig. 4. The neighboring sites of each sublattice are connected by a 
bond if it does not intersect the lines of elementary cell. Each connected cluster of bonds 
in the bulk of lattice is situated inside a loop. Each bulk cluster on the black sublattice is 
surrounded by a connected cluster of bonds on the white sublattice and vice versa. 




FIG. 4: The bijection between bonds and elementary cells. 

The horizontal line cutting the loops cuts bond clusters inside them and produces bound- 
ary clusters surrounded by pieces of loops. Every piece has two points of intersections with 
the horizontal hne and corresponds to the hnk between these points. The clusters of bonds 
corresponding to the loop configuration in Fig. 2 are shown in Fig. 5. 

Since we arc interested in the limiting case of large lattice — )> oo, wc take the infinite 
horizontal line and neglect the left and right boundary conditions. We put the origin of the 
lattice at a selected black site on the horizontal line and consider the boundary clusters of 
bonds containing the origin. The boundary clusters can be classified into four types shown 
schematically in Fig.6. The cluster of type (a) consists of a single vertex coinciding with the 
origin. The type (b) represents clusters containing no points of the horizontal line besides 
the origin. The clusters of type (c) contain the origin and one or more boundary vertices 
left or right from the origin. The clusters of type (d) contain the origin and an arbitrary 
number of boundary vertices left and right from the origin. 

The minimal link probabihty Prob{. . . '-^ . . .) is the sum of probabihties Pa of the clusters 
of type (a) and Pf, of the clusters of type (b) because the minimal link corresponds to the 
piece of loop surrounding the origin (0, 0) together with the cluster of bonds attached to the 
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(-2,0) (0,0) (2,0) 



FIG. 5: The bond configuration corresponding to the loop configuration from Fig. 2. 



(a) 



(b) 





(c) 



(d) 



FIG. 6: Four types of the boundary clusters. 



origin. 



In the 0(1) model, the probability Pa is elementary because the isolated site (0,0) of 
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the black sublattice corresponds to two bonds on the dual white sublattice: one bond con- 
necting sites (—1,0) and (0,1) and another one connecting sites (0,1) and (1,0). For the 
bond percolation, where the bond probability is 1/2, Pa = 1/4. The conjectured value of 
Proh{. . . . . .) is 3/8, then Pb = 1/8. Considering the system of links as a height profile 

[6], one can note that the minimal links correspond to peaks of the profile and the local 
link configurations produced by clusters of type (d) correspond to valleys of the profile. The 
numbers of peaks and valleys in the periodic system coincide, so we have P^ = 3/8. Using 
the identity 

Pa + Pb + 2P, + Pd = l (2.1) 

we conclude that Pc = 1/8. The symmetry P^ = Pc is wonderful because the bond clusters 
of type (b),(c) do not obey any visible symmetry in the percolation theory. In Table 1, 
we show the results of enumeration of clusters of type (a),(b),(c) and (d) for the strip of 
increasing size. We see from Fig. 7 that the convergence of + A to 3/8 is of an order of 
as it is expected from the conjecture (1.1). 
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FIG. 7: Monte-Carlo simulation result. Number of samples is 10*, the N x N upper-left half-plane 
with = 3, 5, 7, . . . , 41 is considered. 

Our aim in this work is calculation of probabilities Pa, Pb, Pc, Pd for the free-fermion 
0(0) model. Elimination of loops in the limit n converts the wavy lines of the 0{n) 
model into dense polymers [9] and the bond clusters on the black and white sublattices into 
spanning trees on these sublattices. A spanning tree on the black (white) sublattice is a 
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N 


Pd 


Pa + Pb 


Pc 


3 


0.0937471400 


0.5313127300 


0.1874740000 


5 


0.1675092500 


0.4423504000 


0.1950355800 


7 


0.2124676800 


0.4119452400 


0.1877856200 


9 


0.2418490300 


0.3982314800 


0.1799565400 


11 


0.2623780800 


0.3909652700 


0.1733454800 


13 


0.2774548300 


0.3866752600 


0.1679902400 


15 


0.2890090100 


0.3838207700 


0.1636097200 


17 


0.2981210300 


0.3818828100 


0.1599813100 


19 


0.3055246300 


0.3805543900 


0.1569469100 


21 


0.3114667000 


0.3797553600 


0.1542892800 


23 


0.3166525700 


0.3788883400 


0.1522516400 


25 


0.3211872800 


0.3781053000 


0.1503894000 


27 


0.3247541400 


0.3778676500 


0.1486765700 


29 


0.3281502900 


0.3774419200 


0.1472263800 


31 


0.3308674800 


0.3771801800 


0.1460154300 


33 


0.3334175200 


0.3769390100 


0.1447843200 


35 


0.3357127000 


0.3766440200 


0.1438112700 


37 


0.3377116700 


0.3765149400 


0.1428866100 


39 


0.3395580300 


0.3764257200 


0.1420234400 


41 


0.3412365600 


0.3762179600 


0.1413121400 


oo 


0.37497 3/8 


0.37499 « 3/8 


0.12506 1/8 



TABLE I: The dense 0{1) loop model: values for = oo are obtained from the extrapolation (see 
Fig.7). 

connected cluster containing all black (white) sites. Due to the cut, the connected spanning 
tree splits into separated components attached to the horizontal line. The classification of 
the components (Fig.6) is the same as in the percolation model. For definiteness, we con- 
sider black components and calculate their statistics via probabilities of white components 
surrounding black clusters. Normally, one selects an arbitrary site of the tree as a root. Then 
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all bonds of the tree become oriented towards the root. It is convenient to choose the root 
of tree on the white sublattice at a site below the horizontal hne that cuts the lattice. Then 
a white cluster having k sites on the cut can be oriented by k different ways because every 
boundary site can be connected with the root by a path along the tree. Evaluation of Pa is 
almost as elementary as that for the 0(1) case. Indeed, we have to find the probabilities of 
three configurations of two oriented bonds on the dual white sublattice: 

la. The bonds from (-1, 0) to (0, 1) and from (0, 1) to (1, 0); 

2a. The bonds from (1, 0) to (0, 1) and from (0, 1) to (-1, 0); 

3a. The bonds from (-1, 0) to (0, 1) and from (1, 0) to (0, 1). 

The probability Pa is the sum of these probabilities, Pa = Pia + P2a + P?,a- 

Evaluation of P^ is a more delicate problem. The cluster of type (b) on the black sublattice 
is surrounded by bonds of a cluster on the white sublattice. The bonds are directed and 
their sequence gives continuous paths on the white sublattice of three different forms: 

lb. The self-avoiding path of an arbitrary length which starts at the site (—1,0), visits 
the site (1, 0), either stops at (1, 0) or continues throughout the white sublattice above the 
cut and ends at a white site on the cut. 

2b. The path starts at the site (1,0), visits the site (—1,0), either stops at (—1,0) or 
continues throughout the white sublattice above the cut and ends at a white site on the cut. 

3b. Two self-avoiding paths start at the sites (—1,0) and (1,0), meet at a white site Sq 
above the cut and then the single self-avoiding path continues from sq up to the end at a 
white site on the cut. 

Denoting the probabilities of these situations by Pib, P26, -P36, we have to find Pf, — Pu + 

Pih + Pzb — Pa- 

While the probabilities constituting are local correlation functions and can be cal- 
culated within a usual approach of the perturbed Laplacian, the probabilities in P^ are 
essentially non-local and their calculation needs a special technique. Before formulating 
these technical tools, we chose more convenient coordinates for our problem. 

Both black and white sublattices are square lattices with the elementary cell \/2 x \/2 
that arc turned by 7r/4 with respect to the basic lattice. We chose the standard orientation 
for the sublattices and put the step of these square lattices to be equal 1. In new coordinates, 
the former cut is the diagonal line crossing the origin. The sites of black sublattice belonging 
to cut have coordinates (0, 0), (1,1), (-1,-1), (2,2), (-2,-2),.... The coordinates of 
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FIG. 8: The upper-left half-plane. Open stars are roots of the spanning trees on the white sublat- 
tice. 

white sublattice on the cut are half- integers (see Fig. 8). To avoid comphcation of notations, 
we shift the origin of the white sublattice to the point (0, 0)' shown in Fig.8. 

After the transformation, the spanning trees on the sublattices belong to the upper-left 
half-space with respect to the diagonal cut. Also we can put the roots of the white sublattice 
at the line parallel to the cut and crossing the point (0, 0)' because all sites of this line can 
be connected with the root of the white spanning tree on the white sublattice. 



III. LAPLACIAN AND GENERALIZED KIRCHHOFF THEOREM 

Consider the graph Q — {V,E) with vertex set V and set of bonds E. The vertices are 
sites of the square lattice and an additional point which is the root : V = {sx,y, {x, y) e 
7j^, \x\ < M,\y\ < N} U {-k}. The bonds of E connect neighboring sites of the lattice. We set 
the condition that all right and bottom boundary vertices {s^^+n} and {s^M,y} are connected 
with the root The graph Q represents the finite square lattice of size (2A^+ 1) x (2M + 1). 
We consider also the left- upper half plane V+ = {sx,y, x < y} with open boundary conditions 
at the diagonal sites Vq = {sx,y, x — y}. 

We construct the spanning tree by using the arrow representation. To this end, we attach 
to each vertex i e an arrow directed along one of bonds (z, i') e E incident to it. 
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Each arrow defines a directed bond {i — > i') and each configuration of arrows defines a 
spanning directed graph with set of bonds {{i i') : i E ^\{*}, i' £ V, e E}. A 

sequence of directed bonds (ii, ^2), (^2, ^3), (^3, ^4), • • • , {im-i,im) is called the path of length 
m from the site ii to the site im- This path forms a loop if im = ii- No arrows are attached 
to the root -k, so that it is a sink of directed paths. Spanning tree is a spanning digraph 
without any loops. Our aim is to construct spanning trees, with a prescribed configuration 
of paths between fixed vertices. These configurations will be investigated by means of the 
determinant expansion of the Laplace matrix [11-13]. 

Let the vertices of the set V, be labeled in arbitrary order from 1 to n = |^\{'^}| = 
(2M + 1)(2A^ + 1). Then Laplacian A of size n x n has the elements: 



Zi if i = j, 

— 1 if i, J are nearest neighbors, (3.1) 
otherwise. 



where Zi is the degree of vertex i G The determinant of A is a sum over all 

permutations a of the set {1, 2, . . . , n}: 

det A = ^ Sgn(c7)Ai,^(i) A2,<,(2) • • • An,a(n) , (3.2) 

where Sn is the symmetric group, sgn(o') = ±1 is the signature of permutation a. In 
general, each permutation a E Sn can be factorized into a composition of disjoint cyclic 
permutations, Ci, C2, . . . c^. This representation partitions the set of vertices into non- 

empty disjoint subsets which are orbits Oi — {i'i,i,i'i,2 ■ ■ ■ ,1^1,1^} C ^ of the corresponding 
cycles Ci, i — 1, . . . ,k, at that U^^^Oi — V\{~k} and X]*Li k — n, where k is the length of 
cycle Ci. Orbits consisting of just one element are the fixed points Sfp{a) of the permutation. 
A cycle q of length |cj| = /j > 2 is called a proper cycle. The proper cycles are of even 
length only, hence, the number of proper cycles p defines the signature of the permutation 
(7, that is sgn((7) = (—1)^- Thus (3.2) can be written as follows: 

n [n/2] p 

detA = n^^ + E(-l)' E U^^i''^i(^^)^<^(^i),'^(^^)---^^-'^0,v. n ^^•'(3-3) 

i=l p=l >7={ci,...,Cp} i=l jeSfp{a) 

where is the A;- fold composition of the cyclic permutation Cj of even length li, Vi & Oi{a), 
so that c^~^{vi) ^ c^{vi) and c''(vj) — Vi. The term HiLi ^« equals to the number of all 
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spanning digraphs having the root Each of other terms on the right-hand side of (3.3) 
having a non-zero set of fixed points equals up to a sign to IljeS/ (a) because all non- 
diagonal elements equal to —1. That product represents the number of distinct spanning 
digraphs which have in common the specified cycles Ci, . . . ,Cp, and differ in the oriented 
edges outgoing from vertices j G Sfp{a). The expansion (3.3) can be interpreted in form 
of the inclusion-exclusion principle [20]. Let Ci, C2, . . . , be the list of all possible proper 
cycles. We define Aj, i — 1,2, ...,m as the set of all spanning digraphs containing the 
particular cycle Cj and Aq is the set of all spanning digraphs. Let ^45^ be the set of spanning 
trees i.e. the set of spanning digraphs containing no cycles. Then we can write down (3.3) 
in the form: 

m 

detA= |^5t| = l^ol Yl \Ar^Aj\ + ... + (-i)^\Ain...nAJ, (3A) 

1=1 l<i<j<m 

where 1^41 is cardinality of the set A. (3.4) is the Kirchhoff theorem for the number of 
spanning trees of a given graph [20] . 

Now we modify the Laplace matrix changing some non-diagonal elements: 

Zi if i = j, 

— 1 if i,j are nearest neighbors, 

(3.5) 

-e if (i, j) e B = {(ii, ji), (i2,j2),---(«n,Jn)}, 



A' = <^ 



otherwise, 

where B is the set of directed bonds that are not nearest neighbors on the square lattice. The 
determinant expansion (3.3) with A' generalizes the Kirchhoff theorem [12]. The product 
A-j . . . A^^ = {~^)^ survives in the limit hm£_^oo det A'/e'^. Permutations corresponding 
to these terms, contain cycles with directed bonds B. If the configuration of bonds B excludes 
cycles containing more than one bond from the set B, the expression lime_^oo det A'/(— e)" 
gives the number of configurations with following features: (i) each configuration is the 
(n + 1)- component spanning graph; (ii) k-th component consists of the path from site 
to site ijt, 1 < /c < n, and branches of the spanning tree attached to these paths; (iii) 
{n + l)-th component is the spanning tree containing the root The ratio of numbers of 
such configurations and all one-component spanning trees is: 

£-^.oo (-£)"det A e-^oo (-s)" 
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where 5 = A' — A is the defect matrix and G — A ^ — G(j,^^q^)^(j,^^q^) is the Green function, 
which is defined in thermodynamical hmit M, A'^ — > oo on the plane as: 

= Go,o+// ^- 3.7 

JJ-n Stt^ 2 — COS a — cos p 

If B includes m bonds {i — > i') between nearest neighbors on the square lattice, the 
diagonal elements A^,j are changed from Zi to e and non-diagonal ones A'.., from —1 to —e. 
In this case, the hmit hm£_).oo det A'/ enumerates the spanning trees with an obligatory 
presence of bonds {i i'). 

IV. CLUSTER PROBABILITIES 

Upon above preparations, we are ready to calculate probabilities Pa — Pia + P2a + -fsa 
and = Pi6 + P26 + Pzh - Pa- 

A. The probability Pa 

The two-bond configuration corresponding to Pa is shown in Fig.9. We denote by B the 
non-zero submatrix of the defect matrix 5. For the case la, matrix B is 

/ 

B{{1 ^ 3,3^ 2}) = 



V 



e 





-A 














—£ 


' ) 



(4.1) 



Pi„ = lim '\ / ' , (4.2) 



The probability Pia follows from (3.6) 

dct(/ + Z?G'"i^) 

where we changed 5 by its non-zero finite submatrix B and the bulk Green function G by 
the Green function (7°^ for the open boundary conditions 

^(Pl,9l),(P2,g2) ~ ^P2-Pl, 92-91 ~ ^P2-9l,g2-Pl ■ (4-3) 

Using exact values (5.6) for the Green functions at short distances given in Appendix, we 
obtain 

167r- 16-37r2 

= 65^ ■ 
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FIG. 9: Two-bond configuration. 



For the case 3a we have 



B{{1 ^ 3,2 ^ 3}) 



^ s -s^ 

£ -e 




(4.5) 



and 



The resulting probabihty Pa is 



3a 



2 (771 - 12 - TT^ 



P. ^ Pa. + 2P.^ii^^^^^ 0.181903, 



(4.6) 



(4.7) 



instead of = 1/4 for the dense 0(1) loop model. 



B. The probability 



The probability Pib is the sum of two terms Pifj{short) and Pib{long). The first of them 
is the probability that the self-avoiding path starting at site (—1,0) of the original white 
sublattice stops at site (1,0). The configuration of this path on the rotated sublattice is 
shown in Fig. 10, where 1 on the original (non-rotated) lattice is the site (—1, 0) and 2 is the 
site (1,0). Two bold arrows denote two possible connections of the site 2 with the roots. 
To generate the path from 1 to 2, wc introduce the bridge (2 — )■ 1) with the weight —e 
(see Fig. 10) and tend e to infinity. Then the defect matrix B consists of the single element 
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FIG. 10: The configuration corresponding to Pih{short) and its bridge representation. 

(2 1) and the formula (3.6) is reduced to 

Pihishort) = 2 hm det(/ + BG^) ^ ^^^^^ _ q^^^^ (4 3) 



Using the table (5.6) we get 



Pu,{short) 



2{n - 3) 



TT 



(4.9) 



The probability Pib{long) involves the self-avoiding paths which start at the site (—1, 0), 
visit (1,0) and then continue through the white sublattice up to the stop at a white site on 
the cut. This situation on the rotated sublattice is depicted in Fig. 11, where points 1, 2, 3, 




it 



FIG. 11: The configuration corresponding to Pib{long) and its bridge representation. 



4 with the original coordinates are correspondingly (—1,0), (1,0), (0, 1), (2, 1). The path 
from the point 4 continues up to its end at a white site on the cut. We generate this path as 
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a sum of two parts: one path from 1 to 2 via 3 and then the path from 4 to the ending point 
5 (see Fig. 11). Each of two parts appears due to auxihary bridges (2 — > 1) and (4 — > 5). 
The orientations of the bridges ensure that a single path through two auxihary bridges is 
forbidden. Then the formula (3.6) can be applied with the result 

/'..(fcn9) = ^^^^^^, (4.10) 

SO that 

Pib = Pibishort) + Pibilong) = ^^"^ ~ ~ . (4.11) 

By the symmetry, we have Pi;, = P2b- 

The crucial point of our calculations is the probability Pst,. The configuration of paths 
on the white sublattice corresponding to the case 36 gets after the transformation the form 
depicted in Fig. 12, where the point 1 is the transformed position of the site (—1,0), the 




FIG. 12: The configuration corresponding to case 36. 

point 2 of the site (1, 0) and the point 4 is the end of a path from the connection point Sq. In 
principle, it is possible to define a system of auxiliary bridges giving all three paths on Fig. 12, 
however, in this case we would obtain a product of three Green functions which should be 
summed up over all positions of the point sq- We choose here another way, avoiding these 
tremendous calculations. 
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First, we define the path from point 1 to point 4 inserting the auxihary bond (4 — > 1). 
Then, we fix an outgoing bond from point 2 in one of two possible directions (diagram I 
shown in Fig. 13). The path from point 2 can join the path from point 1 to 4, forming 




FIG. 13: The diagram I representing Pj. 



thereby a configuration of the type 3b. Otherwise, the path from 2 can finish at one of 
the white sites on the boundary between points 2 and 4. The latter case, that is shown on 
the diagram II, Fig. 14, must be excluded from all configurations in the diagram I. Thus, 
instead of summation over two-dimensional positions of the connection point Sq, we get two 
one-dimensional summations over positions of points 3 and 4, and reduce by one the number 
of Green functions involved into calculations. 

Consider now the Green function G°^^-^ (kk+i) between an arbitrary site {p,q), q > p on 
the left-upper half part of the square lattice and a site {k, k + 1) on the diagonal boundary. 
The summation over the whole boundary gives 



XI ^°p,q),{k,k+l) ~ ^ iGk-p,k+l-q-Gk-q,k+l-p] - (4.12) 

k=—oo k=—oo \ 



Stt^ / / ^ 2 — cos a — COS B 

k = — OD 



r ie sin[(g -p)Qi] ^ J_ T sin a sin[(g - p)Qi] ^ 1 
47r 1 — cos a An J _^ 1 — cos a 2 ' ' 
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FIG. 14: The diagram 11 representing P77. 
With p = and q = 1 we have 

Dk = G'(o,i),(fc,fc+i) = ^'^.fe ~ , (4.14) 

and 

+00 ^ 

E = 2 ■ ^^-^^^ 

A;=— 00 

Due to the symmetry relation 

D.k = Dk. (4.16) 



we can find the sum 

+00 



11 11 
X:B. = j--(G.,„-G„) = j-:^. (4.17) 

fe=l 



Denote by Rk the probabihty of the diagram I for distance k between points 1 and 4. 
The calculation of the determinants with defect matrices B containing the auxiliary bond 
between points 1 and 4, and selected outgoing bonds from the point 2 gives 

+ (4.18) 

TT TT 

The probability Pj is the sum 

00 

Pi = Y,Rk. (4.19) 

k=2 
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Prom the identity 



+00 



Y,R, = 2Y,Rk + ^^^D, + iD_, + Do + D,) (4.20) 



fc = — 00 



k=2 



we have 



67r-8-7r2 



(4.21) 



The probabihty Pu of the diagram II (see Fig. 14) is a double sum of double bridges with 
distance s between points 1 and 4 and distance k between points 1 and 3. The corresponding 
B matrix reads 

/ -e\ 
-e 

\0 / 



B 



which gives the probability 



Pii = 4 hm 

£-^00 



det{I + DG"") 



4 (Dk-iDs - DkDs-i] 



2<k<s 
00 00 



k=2 



2<fc<s fc=2 s=k 

00 00 

= 4 J] L>,(L>i -D,)-AJ2 DkDk-i = 

s=2 ifc=2 
00 00 

= 4Di J] - 4 J] L'ifc(L>fc + Dk-i) , 

s=2 

or, in a more convenient form, 

Pn = 4^1 ^ Dfc - 4 ^ Dk{Dk + Dk+i) . 

k=l k=l 

One can show that 

+00 +00 

J2 Dk{Dk + Dk+i) = 2 J2 MDk + Dk+i) + Do{D^ + 2L>i) 



+00 



k=—oo 



k=l 



SO 



+00 



+00 



Pii = 4Di J] + 2Do{Do + 2D,) - 2 J] Dk{Dk + D^+i) 



k=l 



k=—oo 



Denote the latter sum by Q: 



+00 



(4.22) 



(4.23) 
(4.24) 
(4.25) 
(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



fe=— 00 
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Using the integral representation of the Green function and definition of D^, after some 
manipulations we can show that 

1 ffr (sinai - sin/3i)(sinai - sin(2«i + 2^2 + , ,„ , 
v = t/// 77, ^Tun 7 r- dai d^i dQ;2 . (4.31) 

16 TT'^ JJJ^T, (2 — coscKi — cospi)(2 — cos 0:2 — cos(q;i + pi + 02)) 



By the substitution 



ai Pi a-2 , . 

ti = tanY, t2 = tany, i3 = tany, (4.32) 



sinai = ^^p^, sin^i = ^^^, sina2 = j^^, (4.33) 



cosai = Y^p^, cos/3i = Y^p^, cosa2 = Y^p^, (4.34) 

2 2 2 

dai d/3i dQ;2 = , .2 i , ^2 1 , ^2 "^^2 dta (4.35) 

and symmetrization by permutation t\ i2, we come to the triple integral of a rational 
function 

Q = —All n{ti,h,h)T2{ti,t2,h)dtidt2dh, (4.36) 



^3 



where 



^ . , , , . _ (^1 - ^2)'(tlt2 - l)(tit2 + tits + t2t3 - I)(tit2t3 - ti - ^2 " ^3) 

' - (1 + mi + mi + mti ^m) ' ^ ' 

, 2tit2t3-2t3 + (tl+t2)(t^-l) 

2111, 12, 13) ^ ^ ^_^^2 _ 2tit2t3(tl + t2 + ts) + 2tltltl + (1 + t? + ti)ti ■ ^ > 

After integration over and several manipulations we have 



2 r+oor+t, (ti - t2f {l + tl + tl + tltl - (ti + t2) V(l + t?)(l + ti)) 



^ TT^io (l + t?)^(l + t^)2(t? + ti + 2t?ti) ^^^^^^ 



(4.39) 



The integration over t2 gives 



1 .+00 / 8tf arctan (^1 + 2t?) 4(1 _ ^2^ arctan(ti) 2{n - 2)f\ 



from which we obtain 

4 — TT 

Q = (4.41) 



27r 

As a result, we have 



27r2 - 57r - 4 



77 



TT 



2 



(4.42) 
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and 

P,,^2xPj-2xPjj^ '-^-^ . (4.43) 

TT 

The final result for is 

P, = Pib + P2b + Psb -Pa^ 327r - 8 - Qtt^ ^^^^^ 

The numerical value of this probability is P5 = 1/8 + 0.0001156... that is remarkably close 
to the conjectured value Pb = 1/8 for the dense 0(1) loop model. 

The calculation of the minimal link probability for n = is a first step in a program 
of investigations of more general link correlation functions. Being realized, this program 
would help to establish a correspondence between lattice calculations and operators of the 
logarithmic conformal field theory (LCFT) with the central charge c — —2. Then, the 
operators of the LCFT corresponding to link variables in the dense 0(1) loop model with 
the central charge c = can be constructed by an analogy with the 0(0) model. 

V. APPENDIX: GREEN FUNCTIONS 

The explicit form of the translationally invariant Green function on the plane is 

Ofi,f2 = G{f2 - fi) = Oo,o + gp,q, fa - fi = f = (p, q) (5.1) 
with Go,o an irrelevant infinite constant. The (finite) numbers Qp^q are given explicitly by 

1 rf^ pipa+ig/3 _ i 

9p,, = ^ ^ dad/3. (5.2) 

87r^ 2 — cos a — cos p 

Let us mention symmetry properties of this function: 

9p,q — 9q,p — 9-p,q — 9p,-q ■ (5-3) 

After the integration over a, it can be expressed in a more convenient form for actual 
calculations, 

where t — y — ^/y'^ — 1, y — 2 — cos /3 (see [13]). For = + » 1 it has a behavior [18] 

1 A 3, \ cos(4(/?) 18cos(4(^) + 25cos(8(/?) 
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where {p,q) — (r cos </7,r sine/?), and with 7 = 0.57721... the Euler constant. Let us also 
write down some values of the Green functions for particular points (p, q) , which are needed 
for calculations: 

90.1 — -4 90,2 — -i + - 90,3 — + V 

9i,i^-^ 9i,2^\-l 5'i,3 = 2-|| (5.6) 

92.2 — 92,3 ~ ~4 ~ Stt 93,3 — ~i5^ ■ 
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